The decline of COVID-19 severity and lethality over two years of pandemic

Undernotification of SARS-CoV-2 infections has been a major obstacle to the tracking of critical quantities such as infection attack rates and the probability of severe and lethal outcomes. We use a model of SARS-CoV-2 transmission and vaccination informed by epidemiological and genomic surveillance data to estimate the number of daily infections occurred in Italy in the first two years of pandemic. We estimate that the attack rate of ancestral lineages, Alpha, and Delta were in a similar range (10-17%, range of 95% CI: 7-23%), while that of Omicron until February 20, 2022, was remarkably higher (51%, 95%CI: 33-70%). The combined effect of vaccination, immunity from natural infection, change in variant features, and improved patient management massively reduced the probabilities of hospitalization, admission to intensive care, and death given infection, with 20 to 40-fold reductions during the period of dominance of Omicron compared to the initial acute phase.

Tracking the progressive spread of the SARS-CoV-2 Omicron variant in Italy, December 2021 - January 2022

The SARS-CoV-2 variant of concern Omicron was first detected in Italy in November 2021. Data from three genomic surveys conducted in Italy between December 2021 and January 2022 suggest that Omicron became dominant in less than one month (prevalence on January 3: 78.6%-83.8%) with a doubling time of 2.7-3.1 days. The mean net reproduction number rose from about 1.15 in absence of Omicron to a peak of 1.83 for symptomatic cases and 1.33 for hospitalized cases, while it remained stable for critical cases.

Co-circulation of SARS-CoV-2 variants B.1.1.7 and P.1

SARS-CoV-2 variants of concern (B.1.1.7, P.1 and B.1.351) have emerged in different continents of the world. To date, little information is available on their ecological interactions. Based on two genomic surveillance surveys conducted on February 18 and March 18, 2021 across the whole Italian territory and covering over 3,000 clinical samples, we found significant co-circulation of B.1.1.7 and P.1. We showed that B.1.1.7 was already dominant on February 18 in a majority of regions/autonomous provinces (national prevalence 54%) and almost completely replaced historical lineages by March 18 (dominant in all regions/autonomous provinces, national prevalence 86%). At the same time, we found a substantial proportion of cases of the P.1 lineage on February 18, almost exclusively in Central Italy (with an overall prevalence in the macro-area of 18%), which remained at similar values on March 18, suggesting the inability by this lineage to outcompete B.1.1.7. Only 9 cases from variant B.1.351 were identified in the two surveys. At the national level, we estimated a mean relative transmissibility of B.1.1.7 (compared to historical lineages) ranging between 1.55 and 1.57 (with confidence intervals between 1.45 and 1.66). The relative transmissibility of P.1 estimated at the national level varied according to the assumed degree of cross-protection granted by infection with other lineages and ranged from 1.12 (95%CI 1.03-1.23) in the case of complete immune evasion by P.1 to 1.39 (95%CI 1.26-1.56) in the case of complete cross-protection. These observations may have important consequences on the assessment of future pandemic scenarios.

Return to normal: COVID-19 vaccination under mitigation measures

Vaccination campaigns against COVID-19 are allowing the progressive release of physical distancing restrictions in many countries. However, the global spread of the highly transmissible Delta variant has likely suppressed the residual chances of SARS-CoV-2 elimination through herd immunity alone. Here we assess the impact of the vaccination program in Italy since its start on December 27, 2020 and evaluate possible prospects for reopening the society while at the same time keeping COVID-19 under control. To this aim, we propose a mathematical modeling framework where levels of social activity are adjusted to match the time-series of the net reproduction number as estimated from surveillance data. We compared the estimated level of social contacts, number of deaths, and transmission potential with those of a counterfactual scenario where the same epidemic trajectory is obtained in absence of vaccination. We then evaluate the prospective impact of different scenarios of vaccination coverage and different social activity levels on SARS-CoV-2 reproduction number. We estimate that by June 30, 2021, the COVID-19 vaccination program allowed the resumption of about half the social contacts that were recorded in pre-pandemic times; in absence of vaccination, only about one third could have been resumed to obtain the same number of cases, with the added cost of about 12,100 (95%CI: 6,600-21,000) extra deaths (+27%; 95%CI: 15-47%) between December 27, 2020 and June 30, 2021. We show that the negative effect of the Delta variant diffusion in July was entirely offset by vaccination in the month of July and August 2021. Finally, we estimate that a complete return to the pre-pandemic life could be safely attained only if >90%, including children from 5 years on, will be vaccinated using mRNA vaccines developed in 2020. In any case, increasing the vaccination coverage will allow further margins for societal reopening even in absence of a pediatric vaccine. These results may support the definition of vaccination targets for countries that have already achieved a broad population coverage.

Rapid screening for SARS-CoV-2 variants of concern in clinical and environmental samples using nested RT-PCR assays targeting key mutations of the Spike protein

New SARS-CoV-2 mutations are constantly emerging, raising concerns of increased transmissibility, virulence or escape from host immune response. We describe a nested RT-PCR assay ([~]1500 bps) to detect multiple key spike protein mutations distinctive of the major known circulating SARS-CoV-2 variants, including the three Variants of Concern (VOCs) 20I/501Y.V1 (United Kingdom), 20H/501Y.V2 (South Africa), and 20J/501Y.V3 (Brazil), as well as the 20E.EU1 variant (Spain), the CAL.20C recently identified in California, and the mink-associated variant (GR, lineage B.1.1.298). Prior to application to field samples, the discriminatory potential of this PCR assay was explored using GISAID and Nextclade. To extend variant detection to challenging matrices such as sewage, where the amplification of long fragments is problematic, two short nested RT-PCR assays ([~]300 bps) were also designed, targeting portions of the region spanned by the long nested assay. The three newly-designed assays were then tested on field samples, including 7 fully-sequenced viral isolates from swab samples and 34 urban wastewater samples, some of which collected in areas where circulation of VOCs had been reported. The long assay successfully amplified all the previously characterized viral isolates, allowing the correct identification of variants 20I/501Y.V1 and 20E.EU1 present in the panel. The sequences obtained using the short assays were consistent with those obtained with the long assay. Mutations characteristic of VOCs (UK and Brazilian variant) and of other variant (Spanish) were detected in sewage samples. To our knowledge, this is the first evidence of the presence of sequences harboring key mutations of 20I/501Y.V1 and 20J/501Y.V3 in urban wastewaters, highlighting the potential contribution of wastewater surveillance to explore SARS-CoV-2 diversity. The developed nested RT-PCR assays can be used as an initial rapid screening test to select clinical samples containing mutations of interest. This can speed up diagnosis and optimize resources since it allows full genome sequencing to be done only on clinically relevant specimens. The assays can be also employed for a rapid and cost-effective detection of VOCs or other variants in sewage for the purposes of wastewater-based epidemiology. The approach proposed here can be used to better understand SARS-CoV-2 variant diversity, geographic distribution and impact worldwide.

Effectiveness of regional restrictions in reducing SARS-CoV-2 transmission during the second wave of COVID-19, Italy.

To counter the second COVID-19 wave in autumn 2020, the Italian government introduced a system of physical distancing measures organized in progressively restrictive tiers (coded as yellow, orange, and red) and imposed on a regional basis according to epidemiological risk assessments. The individuals attendance to locations outside the residential settings was progressively reduced with tiers, but less than during the national lockdown against the first COVID-19 wave in the spring. The reproduction number Rt decreased below the epidemic threshold in 85 out of 107 provinces after the introduction of the tier system, reaching average values of about 0.99, 0.89 and 0.77 in the yellow, orange and red tier, respectively. We estimate that the reduced transmissibility resulted in averting about 37% of the hospitalizations between November 5 and November 25, 2020. These results are instrumental to inform public health efforts aimed at preventing future resurgence of cases.

Longevity of seropositivity and neutralizing titers among SARS-CoV-2 infected individuals after 4 months from baseline: a population-based study in the province of Trento

BackgroundThere are conflicting results about the duration of antibodies induced by SARS-CoV-2, but several studies show a rapid decay in a few months after infection. To evaluate antibody decline, we re-evaluated the presence of anti-SARS-CoV-2 antibodies among individuals found seropositive in a first population survey conducted 4 months before. MethodsAll individuals above ten years of age resident in 5 municipalities of the Autonomous Province of Trento, northern Italy, who resulted IgG positive for anti-SARS-CoV-2 nucleocapsid (NC) antibodies in a serosurvey conducted on May 2020 were retested after 4 months. Anti-SARS-CoV-2 antibodies were detected using the Abbott SARS-CoV-2 IgG assay (Abbott Diagnostics, USA) detecting anti-NC antibodies. Samples that gave a negative result were re-tested using the same test plus Liaison SARS-CoV-2 IgG assay (DiaSorin, Italy) to assess anti-spike (S) S1/S2 IgG antibodies. Seroprevalence was calculated as the proportion of positive people on the total number of tested. A neutralizing assay was performed on a subgroup of formerly positives sera using fifty-percent tissue culture infective dose (TCID50) as endpoint dilution to produce a cytopathic effect in 50% of inoculated Vero E6 cells culture. In all the analyses a p value < 0.05 were considered statistically significant. Statistical analysis was performed by STATA version 16.1 (STATA Corp., College Station, Texas, USA). FindingsOverall, 1159 out of 1402 initially anti-NC seropositive participants were enrolled in the study. Of them, 480 (41.1%) became seronegative for anti-NC IgG antibodies. When 479 negative sera were tested for anti-S IgG, 373 samples (77.9%) resulted positives. A functional neutralization assay was performed on 106 sera showing high concordance with anti-S antibodies positivity. InterpretationA decline of anti-NC IgG values was recorded 4 months after the first evaluation. Worth of note, a high proportion of anti-NC seronegative individuals were positive for anti-spike IgG antibodies, which appear to persist longer and to better correlate with neutralization activity.

Epidemiological characteristics of COVID-19 cases in non-Italian nationals notified to the Italian surveillance system

BackgroundInternational literature suggests that disadvantaged groups are at higher risk of morbidity and mortality from SARS-CoV-2 infection due to poorer living/working conditions and barriers to healthcare access. Yet, to date, there is no evidence of this disproportionate impact on non-national individuals, including economic migrants, short-term travellers, and refugees. MethodsWe analysed data from the Italian surveillance system of all COVID-19 laboratory-confirmed cases tested positive from the beginning of the outbreak (20th of February) to the 19th of July 2020. We used multilevel negative-binomial regression models to compare the case-fatality rate and the rate of admission to hospital and intensive care unit (ICU) between Italian and non-Italian nationals. The analysis was adjusted for differences in demographic characteristics, pre-existing comorbidities, and period of diagnosis. ResultsWe analysed 213,180 COVID-19 cases, including 15,974 (7.5%) non-Italian nationals. We found that, compared to Italian cases, non-Italian cases were diagnosed at a later date and were more likely to be hospitalised [(adjusted relative risk (ARR)=1.39, 95% confidence interval (CI): 1.33-1.44)] and admitted to ICU (ARR=1.19, 95% CI: 1.07-1.32), with differences being more pronounced in those coming from countries with lower HDI. We also observed an increased risk of death in non-Italian cases from low-HDI countries (ARR=1.32, 95% CI: 1.01-1.75). ConclusionsA delayed diagnosis in non-Italian cases could explain their worse outcomes compared to Italian cases. Ensuring early access to diagnosis and treatment to non-Italians could facilitate the control of SARS-CoV-2 transmission and improve health outcomes in all people living in Italy, regardless of nationality.

Epidemiological characteristics of COVID-19 cases in Italy and estimates of the reproductive numbers one month into the epidemic

BackgroundIn February 2020, a locally-acquired COVID-19 case was detected in Lombardia, Italy. This was the first signal of ongoing transmission of SARS-CoV-2 in the country. The outbreak rapidly escalated to a national level epidemic, amid the WHO declaration of a pandemic. MethodsWe analysed data from the national case-based integrated surveillance system of all RT-PCR confirmed COVID-19 infections as of March 24th 2020, collected from all Italian regions and autonomous provinces. Here we provide a descriptive epidemiological summary on the first 62,843 COVID-19 cases in Italy as well as estimates of the basic and net reproductive numbers by region. FindingsOf the 62,843 cases of COVID-19 analysed, 71.6% were reported from three Regions (Lombardia, Veneto and Emilia-Romagna). All cases reported after February 20th were locally acquired. Estimates of R0 varied between 2.5 (95%CI: 2.18-2.83) in Toscana and 3 (95%CI: 2.68-3.33) in Lazio, with epidemic doubling time of 3.2 days (95%CI: 2.3-5.2) and 2.9 days (95%CI: 2.2-4.3), respectively. The net reproduction number showed a decreasing trend starting around February 20-25, 2020 in Northern regions. Notably, 5,760 cases were reported among health care workers. Of the 5,541 reported COVID-19 associated deaths, 49% occurred in people aged 80 years or above with an overall crude CFR of 8.8%. Male sex and age were independent risk factors for COVID-19 death. InterpretationThe COVID-19 infection in Italy emerged with a clustering onset similar to the one described in Wuhan, China and likewise showed worse outcomes in older males with comorbidities. Initial R0 at 2.96 in Lombardia, explains the high case-load and rapid geographical spread observed. Overall Rt in Italian regions is currently decreasing albeit with large diversities across the country, supporting the importance of combined non-pharmacological control measures. Fundingroutine institutional funding was used to perform this work.

Appropriate laboratory utilization in diagnosing pulmonary embolism.

BACKGROUND: The appropriateness of clinical laboratory use in hospital clinical practice is a very debated question. In fact diagnostic tests are essential tools for disease screening or diagnosis but at the same time they represent an important expenditure. METHODS: A quantitative study was conducted in an Italian general hospital with 535 beds and about 27,000 admissions per year. The sample was made of all patients discharged from the hospital with DRG 78 between the period 1 January 2005 and 31 December 2005. RESULTS: The Emergency Department (ED) discharged 2.9% (116/4009) of patients with pulmonary embolism diagnosis in the year 2005. The percentage of prescription inappropriateness inferred by analysis of all required tests by operative unit was 21.7% (950/4385). The approximate estimate of the economic value of unnecessarily required tests was of 3495 euro. Of these 96.5% (112/116) had enlisting criteria. All haematological and clinical-chemical tests (1295) concerning studied patients were analysed. About 70.4% (93/132) of d-dimer tests were ordered in ED. In the studied patients, unfractionated heparin was administered in 17.8% (20/112) of the cases, low-molecular weight heparin in 79.5% (89/112) and heparin therapy was not administered in 2.7% (three of 112). CONCLUSION: This study uses a method to assess the quality of laboratory test orders using results as a tool to estimate the impact of economic resources devolved in executing inappropriate tests.